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ABSTRACT 



Aims. To calculate chemistry and gas temperature of evolving protoplanetary disks with decreasing mass or dust settling, and to 
explore the sensitivity of gas-phase tracers. 

Methods. The density and dust temperature profiles for a range of models of flaring and self-shadowed disks around a typical Herbig 
Ae star are used together with 2-dimensional ultraviolet (UV) radiative transfer to calculate the chemistry and gas temperature. In 
each model the line profiles and intensities for the fine structure lines of [O i], [C n] and [C i] and the pure rotational lines of CO, CN, 
HCN and HCO + are determined. 

Results. The chemistry shows a strong correlation with disk mass. Molecules that are easily dissociated, like HCN, require high 
densities and large extinctions before they can become abundant. The products of photodissociation, like CN and C2H, become 
abundant in models with lower masses. Dust settling mainly affects the gas temperature, and thus high temperature tracers like the 
O and C + fine structure lines. The carbon chemistry is found to be very sensitive to the adopted PAH abundance. The line ratios 
CO/ 13 CO, CO/HCO + and [Oi] 63yum/146yum can be used to distinguish between disks where dust growth and settling takes place, 
and disks that undergo overall mass loss. 
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1. Introduction 

Circumstellar disks are a natural by-product of the formation 
of stars of low to intermediate masses, the former objects be- 
ing identified as T-Tauri stars while the latter are Herbig Ae/Be 
stars. From models of the spectral energy distributions (SEDs) at 
infrared and millimeter wavelengths, which track the dust emis- 
sion from these disks, it has been concluded that some disks have 
a flaring geometry, allowing them to intercept a significant por- 
tion of the st ellar radiation and re-emitting it at longer wave- 
length s (e.g. iKenvon & Hartmannl I1987I: IChiang & Goldreichl 
Il997l) . In a recent study, iDullemond & Dominikl (|2004) sug- 
gested that the difference between flaring ('group F) and flat 
('group IF) disks can be explained by the presence of a puffed- 
up inner rim that casts a shadow over the disk. Depending on 
whether there is enough material in the outer disk to rise above 
the shadow either an illuminated flared disk or a self-shadowed 
flat disk results. One interpretation of the occurrence of group 
I and group II sources is that flared disks evolve into flat disks 
when the dust grows and settles to the midplane, or when the 
disk's mass decreases due to dissipation and evaporation. 

While SEDs are useful to examine the properties of the 
dust content of disks, the physical and chemical properties 
of the gas are studied through line emission from both fine- 
structure transitions in atoms and rotational and vibrational tran- 



lines of CO (e.e. iKoerner & Sargent 


1995: Dutreyetal.il 19971 


iThi et all I2001 : Brittain et al.l I20031 


Dent et all 120051). and H 7 



as the rotational lines of several other molecules like CN , 
HCN, HCO + and C 2 H (e.g. lDutrev et~aDfT997l: IThi et al.ll2004b . 
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To obtain information on the temperature, chemical compo- 
sition and total mass of the gas from the observed emission 
lines, the chemistry of the disk has to be mo delled. In re- 
cent years calculations of both th e chemistry (e.g. lAikawa et al.l 
1997; van Zadelhoff et alj|2003l) and the gas te mperature (e.g. 
iJonkheid et al.ll2004t iKamp & Dullemondll2004l) have been per- 
formed using a density structure based on previous models of 
the d isk structure (e.g. iD'Alessio et al.|[l998l: IDullemond et al.l 
2002). In these models, the density structure is calculated as- 
suming vertical hydrostatic equilibrium, and the dust temper- 
ature is obtained by solving the radiative transfer of stel- 
lar light. In most cases the gas temperature is assumed to 
be equal to the dust temperature, but some models calculate 
the gas temperat ure and density structure i n equilibrium with 
each other (e.g. iGorti & Hollenbachl 12004 iNomura & Millar! 
l2005h . The emission lines are modelled from the results of 
the chemistry and temper a ture calculations (e.g. IKamp et al.l 
120031: IJonkheid et alJl200l IGorti & Hollenbachl 120041) . Effects 
of grain growth and dust settling on the chemistry have also been 
considered (e.g. IKamp et all 120031: IJo nkheid et alj|2004 120061: 
lAikawa & Nomura 200^~ 



In the current paper calculations of the chemistry and gas 
temperature are presented based on the disk structure models of 
group I and group II sources by Dullemond & Dominik (2004) 
to investigate how the chemistry and line emisison change in 
a series of evolving disks. Special attention was given to iden- 
tify possible tracers of the disk geometry that can supplement 
and confirm the information gained from SEDs. This work com- 
pleme nts the calculations for optically thin disks by Kam p et all 
(20061). 



Jonkheid et al.: Chemistry and line emission from evolving Herbig Ae disks 



M 100 



SO 




n (em 3 ) 








100 200 300 400 500 
R (AU) 



100 800 300 400 S00 
R (AU) 



100 300 300 400 500 
R (AU) 



Fig. 1. Input densiti es (left column) dust tempera tures (middle column) and integrated strength of the UV field (right column) for 
models B1-B4 of iDullemond & Dominild (|2004|) . The gas mass (in M Q ) and dust/gas ratio are indica ted in the for each row. The 
UV field is given in terms of /uv times the strength of the interstellar radiation field bv lDrainel (1 19781) integrated between 912 and 
2050 A. The dotted lines in the density plots give the boundary where the visual extinction towards the central star is 1 magnitude. 



2. Model 

Th e chemistry is calculated for a r ange of disk structure models 
by |Dullemond & Dominik (2004) in which the disk mass and 
dust/gas ratio are varied systematically (see TableQ]and Figures 
Q]and[2|i. In that paper the vertical structure of the disk is com- 
puted in hydrostatic equilibrium assuming several radial distri- 
butions and total disk masses. In all models the gas temperature 
was taken equal to the dust temperature. Under consideration 
here are models B 1-B4 and BL1-BL4. These models have a sur- 
face density distribution that falls off as S oc R 15 for R < 200 
AU; beyond 200 AU the surface density drops off as S oc R~ 12 
to simulate an effective outer edge of 200 AU while keeping a 
continuous mass distribution. In the B series of models the total 



disk mass is varied from 10~' M in Bl to 10~ 4 M in B4, while 
keeping a constant dust/gas mass ratio of 0.01. In the BL series 
the gas mass stays constant at 10~' M and the dust/gas ratio is 
varied from 0.01 in BL1 to 10" 5 in BL4. IDullemond & Dominild 
(|2004|) kept the overall dust/gas ratio in the BL series constant 
at 0.01 by removing the interstellar-sized grains (~ 0.1 fxm) from 
the disk and replacing them with cm-sized grains in the mid- 
plane, thus simulating the effects of dust settling and growth, 
but these large grains are excluded from the chemical calcula- 
tions presented here because of their very limited i nfluence (e.g. 
iKamp & van ZadelhoflfeOOltUonkheid et al.ll2006l) . It should be 
noted that models Bl and BL1 are identical. The abundance of 
polycyclic aromatic hydrocarbons (PAHs) is assumed to follow 
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Fig. 2. As FigureG] but for models BL1-BL4. 
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Table 1. Model parameters used in the calculations. Note that 
the dust masses given here include only small (interstellar sized) 
grains. 



Table 2. Adopted gas-phase elemental abundances with respect 
to total hydrogen. 



Model 


M gas (M ) 


M dust (M ) 


Bl 


io-' 


1CT J 


B2 


io- 2 


1(T 4 


B3 


1(T 3 


IO" 5 


B4 


IO" 4 


10~ 6 


BL1 


io- 1 


io- J 


BL2 


io- 1 


IO" 4 


BL3 


io- 1 


IO" 5 


BL4 


lcr 1 


lCr 6 



the dust/gas ratio (i.e. the PAH abundance drops off with a fac- 
tor 10 in each step of the BL series), consistent with the find- 



Element 


abundance 


C 


1.2 x 1(T 4 


O 


1.6 x 1(T 4 


N 


2.0 x IO -5 


Mg 


4.2 x 1(T 6 


S 


1.9 x 1(T 6 


Si 


8.0 x 10~ 7 


Fe 


4.3 x 10~ 6 


PAH 


1 x 10-'° - 1 x 10~ 7 a 



In the B series of models, the PAH abundance is kept constant at 
10~ 7 . In the BL series, the PAH abundance is lowered together with 
the dust content of the disk, going from 10~ 7 in BL1 to 10~ 10 in 
BL4. 
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Fig. 3. The radiation fields used in the model. The solid line 
gi ves the s t ellar radiation field, scaled to the interstellar UV field 
by lDrainel (Il978h (dotted line). The horizontal lines indicate the 
wavelength ranges where several important molecules are disso- 
ciated. 



ings of lAcke & van den Anckeri d2004 and lHabart et al.1 d2004 
that PAH emission is weaker in group II sources than in group 
I sources. The abundances of the elements used in the calcula- 
tions, given in Table|2] are representative for disks where N and 
O are mostly depleted in some form (e.g. H2O ice) onto dust 
rains. 

The densities and dust t e mpera tures are taken directly from 
the iDullemond & Dominikl (2004) models. Only regions with 
densities «h > 10 3 cm~ 3 were taken into account, regions with 
lower densities do not contribute enough column density to be 
taken into consideration. Note that the dust temperature does not 
fall below 30 K anywhere in the disk, so the freeze-out of CO can 
be disregarded. The maximum outer radius of the computational 
grid is taken to be 500 AU. 

In the B series (Figure Q]) the disk is initially very large and 
flared. As the mass decreases in this series, the disk becomes 
smaller in both the radial and the vertical direction. Although 
the disk is still rather thick in models B3 and B4, the t = 1 
surface for stellar radiation shows that there is little flaring: the 
disk does not intercept much radiation at large radii. These mod- 
els are a transition between group I sources (of which Bl and 
B2 are examples) and group II sources (B5 and B6), which have 
too low gas masses to be considered in this paper. The structure 
of the dust temperature and the radiation field show that the in- 
ner part of the disk is strongly self-shadowed. As the disk mass 
decreases, this effect becomes less strong and the low dust tem- 
peratures and UV intensities are confined to a narrow region near 
the midplane. 

The BL series is very similar to the B series regarding the 
temperature structure and UV intensities. While one would ex- 
pect more flaring in models BL3 and BL4 (since the gas distri- 
bution is higher in the outer disk), the r = 1 surface again shows 
that these disks have a largely flat geometry. As in the B series, 
these models form a transition between group I (BL1 and BL2) 
and group II (BL5 and BL6) sources. 

The transfer of UV radiation is per f ormed with the 2-D 
Monte Carlo code by Ivan Zadelhoff et al.l ([2003 ) to account for 
the effects of forward scattering and to include the interstellar 



radiation field. The self-shielding of H2 and CO are included 
by assuming constant abundances for these species of 0.5 and 
10~ 5 , respectively; the results do not depend strongly on the val- 
ues used. The shielding factors are calculated with the mean 
column den sity computed by the Mo nte Carlo code (for more 
details, see Ivan Zadelhoff et all 120031) . The results of the UV 
transfer are used to calculate the photorates dir ectly from the 
cross -sections and local radiation field (following lvan Dishoeckl 
1988) and photoelectric heating rates on classical grains and 
PAHs dTielens & Hollenbach|[l98l iBakes & Tielenslll994 re- 
spectively). 

For the stellar radiation field a NE XTGEN spectrum was 
adopted with T e $ = 10000 K, logg = 4 dHauschildt etalJI 19991 
see Figure O and a stella r radius of The interstellar radia- 
tion field was t aken from Draine ( 1978), with the extension for 
A > 2000 A bv Ivan Dishoeck & Blackl d!982h . Since the formu- 
lation of the photoelectric heating rates does not take the detailed 
spectrum of the local radiation field into account, the correction 
factor by Spaans et al. (1994) (Equation 13 in that paper) was 
used to account for the 10 000 K stellar spectrum. The value /uv 
characterizing the strength of the UV field was obtained by inte- 
grating the radiation field betwe en 912 and 2 050 A and dividing 
this by the same integral for the lDrainel d!978l) field. 

With the photorates and photoelectric heating rates known, 
th e chemistry and gas temperature are calculated as outlined 
in iJonkheid et al.l (|2004): the chemical network of iJansen et al.l 
d 19951) is used to calculate the gas-phase chemistry in steady 
state, and the thermal balance is solved using an escape probabil- 
ity formalism assuming cooling radiation escapes from the disk 
in the vertical direction (i.e., column densities used in the escape 
probability are calculated vertically). The heating processes in- 
cluded in the calculations are: photoelectric effect on both clas- 
sical grains and PAHs, cosmic ray ionization of H and H2, H2 
formation, H2 photodissociation, collisional de-excitation of the 
excited vibrational levels of H2 and the fine structure levels of O, 
and C photoionization. The gas is cooled by line emission from 
the fine structure lines of C + , C and O, the rotational lines of CO, 
and thermal coupling with dust grains. A more thorough descri p- 
tion of these processes can be found in lJonkheid et alj d2004l) . 

The line profiles are calculated from the chemistry and 
gas temperature using the 2-D cyl indrical symmetric acceler- 
ated M onte Carlo code RATRAN bv lHogerheijde & van der Takl 
(2000). To speed up the calculations only the material outside 
R = 10 AU is considered, since the inner part of the disk con- 
tributes little to the overall emission of the atomic and molecular 
lines under consideration in a typical submm observing beam 
(~ 15")- The number of grid cells varies between 900 and 1500, 
depending on the model. 



3. Results and discussion 

3.1. Gas temperature 

The calculated gas temperatures are illustrated in Figures [4] and 
[6] for models B1-B4, and in Figures and [7] for models BL1- 
BL4. In the B series, the surface of the disk is very hot (sev- 
eral thousend K) due to the photoelectric heating, which is 
very efficient compared to the cooling by [Oi] and Lye. The 
temperatures found here are consistent with PDR models with 
«h = 10 3 cirT 3 and I^y = 10 4 - 10 5 . The higher densities in 
the optically thin part of the disk below the immediate surface 
increase the efficiency of [0 1] cooling, resulting in lower tem- 
peratures (several hundred K). Deeper into the disk the densities 
are so high (10 7 - 10 s cm 4 ) that gas-dust coupling dominates 
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Fig. 4. Gas temperatures (left column), CO abundances (middle column) and HCN abundances (right column) for models B1-B4. 



the thermal balance and the gas temperature is equal to the dust 
temperature. In the outer parts of the disk (R > 200 AU) the 
densities are not high enough for gas-dust coupling to dominate 
due to the rapid decline of the total column density. In those re- 
gions, the gas temperature drops below the dust temperature due 
to efficient CO cooling. 

It should be noted that th e gas temperatures fou nd here are 
higher than those reported in Jonkhei d et al] (120041) for T-Tauri 
disks. This is because the radiation fields used in the current pa- 
per are almost two orders of magnitude larger (see Figure [TJ: at 
a radius of 100 AU the current models have /uv = 5 - 6 x 10 4 
at the disk's surface, in contrast with a value of 5 - 6 x 10 2 in 
Jonk heid et ajj (2004). T his is due to the m ore luminous central 
star used by Dullemond & Dominik (2004) and this work. 



In the BL series the importance of photoelectric heating de- 
creases with each step because both small dust grains and PAHs 
disappear from the disk. This means that the total heating rate 
is less than in the B series, and consequently the temperatures 
at the disk's surface are lower (20-200 K going from BL4 to 
BL2). In models BL3 and BL4 heating by C ionization becomes 
important, and this heating rate has a similar magnitude as the 
photoelectric processes. At large depths, gas/dust coupling dom- 
inates the thermal balance, even with the low dust/gas ratio in 
these models. Since the gas temperature is lower than the dust 
temperature in much of the disk, this means that the temperature 
increases toward the midplane. 



3.2. Chemistry 

The vertical distribution of chemical species resembles the struc- 
ture commonly found in PDRs: a thin layer of mostly atomic 
hydrogen at the surface, while slightly deeper in the disk H2 be- 
comes self-shielding and thus hydrogen becomes mainly molec- 
ular. At the surface carbon is mostly in the form of C + and con- 
verts to C and CO deeper in the disk. At very high optical depths 
atomic oxygen is locked into O2. 

There are several deviations from the classical PDR struc- 
ture. In the B series there is a layer of relatively high CO abun- 
dance in the optically thin part of the disk (at z > 50 AU), where 
CO is not self-shielding. This is caused by the high tempera- 
tures in this layer, which help formation of CO via the reac- 
tion C + + H 2 -» CH + + H, which has an enery barrier of 4600 K. 
CH + in turn reacts with O to form CO + , which leads to HCO + 
and CO. Although CO is completely self-shielded in the mid- 
plane of most of the models (an exception being the outer re- 
gions in model B4, where the CO photodissociation by the inter- 
stellar UV field is unshielded), the abundance of C is still rela- 
tively high there. This is caused by the shape of the stellar radia- 
tion field (see Figure|3}: there are still many continuum photons 
with A > 1200 A which can dissociate precursor molecules like 
CH, which is an important pathway of CO formation. In the later 
steps of the BL series C is the dominant form of carbon in the 
midplane instead of CO. This is partly due to the photodissoci- 
ation of precursor molecules, and partly because the formation 
of CO is further suppressed due to lack of PAHs in these mod- 



6 



Jonkheid et al.: Chemistry and line emission from evolving Herbig Ae disks 




D 100 200 300 400 900 100 £00 300 400 500 100 ZOO 300 400 500 



sso r 

£00 [ 
150J 

ioo[ 

50 [ 





350 




£00 




ISO 


< 




N 


100 




M 










£50 




£00 




150 


< 




100 




30 








100 £00 300 400 500 



100 ZOO 300 400 800 



3000 




250 


1000 




£00 


300 




ISO 


100 

30 

10 




100 
50 
Q 

( 


3000 




250 


1000 




£00 


300 


-3 


ISO 


100 

30 

10 


S 


100 

SO 


e 



I 



3000 

1000 

300 

100 

30 

10 




100 £00 300 400 500 



100 £00 300 400 500 



250 
£00 
ISO 
100 

so 






£00 



100 
SO 




250 
£00 
ISO 
100 
SO 




100 £00 300 400 500 



100 £00 300 400 500 



1 




250 






10" a 
10-1 


Da 

& 
s 

i 


£00 
150 


10" s 
10-' 


ft 


d 


100 


10 , 


o 


SO 


10"* 








1 10 " 

io- 
\&* 



I 10" 
10" 

io- 

■if 

Iff 



3- 

c 

3 




200 300 400 
R <AU) 



100 200 300 400 500 
H (AU) 



100 £00 300 400 500 
R (AU) 



Fig. 5. Gas temperatures (left column), CO abundances (middle column) and HCN abundances (right column) for models BL1-BL4. 



els. The reaction PAH : H + C + -> PAH + CH + (where PAH:H 
is a PAH with an extra H atom) is another important step in 
the reaction pathw ay leading to CO (for rate coefficient, see 
ICanosa et aDll995h . In the later steps of the BL series the PAH 
abundance is very low, and so this pathway is closed. 

The global distributions of CO and HCN are illustrated in 
Figure [4] for models B1-B4, and in Figure [5] for models BL1- 
BL4. The abundances of CO and HCN decrease significantly 
going from model B 1 to B4, while this effect is less pronounced 
going from BL1 to BL4. Although the UV continuum flux is 
comparable in both series, the total gas mass stays at 0. 1 M in 
the BL series while it decreases with factors of 10 in the B series. 
The self-shielding of CO will therefore be much stronger in the 
BL series, resulting in higher CO abundances. HCN is rapidly 
dissociated in the optically thin regions of both series, but its 
formation rate correlates with CO abundance. Therefore, a rapid 
decrease in HCN abundance can be seen in the B series, where 
there is not enough CO to counter the rapid photodissociation. 
In the BL series, however, the high CO abundances increase the 
production rate of HCN, keeping it at a higher abundance than 
in the B series. 

A more detailed view of the chemistry is given in Figures 
[6] and [7] which show a vertical cut at 200 AU. In the B series 
the high CO abundances at the surface result in high abundances 
of HCO + , C 2 H, CN and HCN. The abundance of C 2 H follows 
the same trends as CN, wh ich is in agreement with the results by 
Ivan Z adelhoff et al. (20031) forT-Tauri stars. Note that in contrast 



withfvan Zadel hoff et alj d2003f) (who found similar abundances 
for these molecules) the abundance of C 2 H found here is one to 
two orders of magnitude less than CN. This is due to the differ- 
ence in stellar spectrum used in the current paper, which has a 
jump in intensity at wavelengths longer than 1200 A. Since C 2 H 
can be dissociated by photons with A > 1200 A and CN is dis- 
sociated only at A < 1 150 A, the photodissociation rate of C 2 H 
is much higher, and thus this molecule has a lower abundance. 

The same effect is visible in the abundance of HCN com- 
pared to CN: HCN is easily dissociated by photons with A > 
1200 A, giving it a lower abundance than CN in most of the 
disk. Only in the heavily shadowed regions near the midplane 
is HCN more abundant than CN. In the BL series the densities 
near the midplane are high enough to increase the HCN forma- 
tion rate, and so the HCN abundance is higher than CN even in 
the optically thin models BL3 and BL4. 

The abundances of HCO + in the deeper layers of the disk 
are lower by several orde rs of magnitude than those found by 
Ivan Zadelhoff et all d2003l) . A possible explanation for this may 
be that charge exchange with metals plays a more important 
role in the current model, since the densities are a factor ~ 100 
higher. Also, as for CO, the rapid photodissociation of hydrocar- 
bon precursors plays a role. 

In summary, the B series shows a clear trend in molecular 
abundances: going from model Bl to B4 molecules are more 
easily dissociated, and the peaks in their vertical distributions oc- 
cur closer to the midplane. In the BL series the trend is less clear. 
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Fig. 6. Vertical profiles at a radius of 200 AU of the gas temperature and density (left column) and chemistry (middle and right 
columns) for models B1-B4. The visual extinction toward the central star is given at the top of each panel. 



Going from model BL1 to BL2, the drop in gas temperature 
greatly affects the carbon chemistry. The differences between 
models BL2 to BL3 are limited, and models BL3 and BL4 are 
nearly identical. Again there is a shift toward the midplane in the 
peak of the vertical abu ndance distribu tion, in agreement with 
earlier findings of Jonkh eid et al.l {2004) and lAikawa & Nomura! 
(120061) . but this shift is small compared to the B series. 



3.2.1. Influence of PAHs 



form PAH + X + — > PAH + + X) and hydrogenation reactions (of 
the form PAH : H + X -> PAH + XH). Especially the latter reac- 
tions are important in the formation of molecules, as can be seen 
in the distribution of CO in the BL model series (seen in Figures 
[5] and |7), where the PAH abundance is decreased in each step. 
As a result, the reaction C + + PAH : H — > CH + + PAH becomes 
less efficient; since CH + is a chemical precursor of CO (via CO + 
and HCO + ), the overall CO abundance decreases together with 
the PAH abundance. 



It is assumed that PAHs are present at interstellar abundances 
(«pah/«h = 10~ 7 ) for a dust/gas ratio of 10~ 2 , and that these 
abundances scale with the dust/ gas ratio. This assumption is 
con sistent with the find ings of lAcke & van den Anckerl ([2004) 
and lHabart et al.l ([2004), although models by Geers et al. (in 
press) suggest an abundance that is lower by as much as an 
order of magnitude. The influence of PAHs on the chemistry 
can best be seen in the BL series (Figure [7]), where the PAH 
abundance is decreased by a factor 10 at each step. PAHs affect 
the chemistry mainly through charge exchange reactions (of the 



3.2.2. Carbon chemistry 

In some of the models, particularly BL3 and BL4, there is a 
broad layer near the disk midplane where atomic carbon is the 
main carbon-bearing species. This somewhat unconventional re- 
sult is caused by shape of the stellar radiation field, shown 
in Figure [3] the stellar radiation has few photons in the 912- 
1200 A range, most of the flux is at longer wavelengths. This 
means that the photoionization rate of C and the photodissoci- 
ation rate of CO are small, since these processes occur in the 
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912-1 100 Arange. The photodissociation rate of CH is relatively 
high, however, since this molecule is dissociated at longer wave- 
lengths than CO. Since CH is an important chemical precur- 
sor to CO, the rapid dissociation of CH means that the forma- 
tion rate of CO is very low, resulting in a low CO abundance. 
Since atomic carbon is released by the photodissociation of most 
carbon-bearing molecules, and the ionization rate of atomic car- 
bon is very low, carbon is mostly neutral and atomic. 

3.3. Column densities 

The vertical column densities of CO, CN, HCN, C 2 H and HCO + 
are shown in Figure [8] It can be seen that the greatest variations 
occur in the B series for CO and HCN. The CO abundance is 
sensitive to the gas mass due to self-shielding. Similarly, HCN 
occurs only in regions that are dense and have sufficient con- 
tinuum optical depth toward the central star. Both of these con- 
ditions are strongly dependent on the disk mass. The column 
densities of CN and C 2 H in the B series show a peak which 
starts at relatively large radii and shifts inward going from B 1 
to B4. These molecules require a moderate optical depth toward 



the central star: large enough so they and their precursors are not 
dissociated, but small enough that they can still be formed by 
photodissociation of less stable molecules (HCN for CN; C2H2 
and C3H for C2H). The HCO + peak also shifts inward due to the 
lower abundance of CO in the outer regions of the disk. 

In the BL series, the CO column densities show a clear dis- 
tinction between on the one side models BL1 and BL2, and on 
the other side BL3 and BL4. The reason for this can be seen in 
Figure |7J in models BL1 and BL2, CO is the main form of car- 
bon in the dense regions near the midplane. In BL3 and BL4, 
C becomes the dominant form of carbon due to the low PAH 
abundance in these models. In the BL series HCN occurs mainly 
in the midplane of the disk, in a layer that gets narrower going 
from BL1 to BL4. This explains the decreasing column densi- 
ties in Figure[8] The peaks in the vertical distributions of CN and 
C2H get broader going from BL1 to BL4, and consequently their 
column densities in the main portion of the disk (R < 200 AU) 
increase. In the outer disk the interstellar radiation field can dis- 
sociate these molecules, decreasing their column densities. The 
peak in the vertical HCO + distribution decreases in height go- 
ing from model BL1 to BL4, which is reflected in the column 
densities. 
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Fig. 8. The vertical column densities of CO (top row), HCN (second row), CN (third row), C2H (fourth row) and HCO + (bottom 
row) for the B series (with a descending disk mass, left column) and the BL series (with a descending dust/gas ratio, right column). 
The solid lines give the results for model B1/BL1; the dotted line for B2/BL2; the dashed line for B3/BL3; the dash-dotted line for 
B4/BL4. 



The column densities of CN and C2H (and to a lesser extent 
those of HCO + ) are relatively high in the outer parts of the disk 
(R > 200 AU). In spite of the fact that the total column density 
drops off as R~ 12 , the densities are still high enough to provide 
high abundances for these molecules. The abundances of these 
molecules in the outer disk depend strongly on the exact form 
of the total column density profile. If it drops off more rapidly 
than R 12 , this will be reflected in the intensities of the molecu- 
lar lines. As a consequence, the emission lines of CN, C2H and 
HCO + may be used to examine the outer regions of disks that 
are difficult to probe in the continuum. 

The trends in the column densities of CO, CN and HCO + 
presented here for the BL series, in which dust growth and set- 



tling is simulated, show q ualitative similarities to those found by 
lAikawa & Nomura (120061) for T-Tauri disks with different stages 
of grain growth. With increasing grain size, they find that ^V(CO) 
changes only very little, N(CN) stays constant in the outer re- 
gions and increases in the inner regions, and A^(HCO + ) decreases 
in the outer regions. They find no significant variation for HCN 
and C2H, whereas column densities of these molecules change 
significantly in the current paper. This is possibly due to the jump 
in the stellar radiation fie ld at 1200 A, which has a strong influ- 
ence on these molecules; I Aikawa & No mura (2006|) use a much 
smoother stellar spectrum that resembles the interstellar radia- 
tion field. 
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Fig. 9. Emission line profiles for O (top left), CO (top right), HCN (lower left) and CN (lower right) for models B 1 /BL 1 (solid line), 
B4 (dotted line) and BL4 (dashed line). 



3.4. Atomic and molecular lines 

From the molecular abundances and the gas temperature at 
each position in the disk the profiles of rotational molecular 
lines are calculated using the acce lerated Monte Carlo code by 
Hogerheii de & van der Takl (|2000). The distance to the disk is 
taken to be 150 pc, and an inclination of 45° is assumed. For 
the velocity field Keplerian rotation around a 2.5 M star and a 
turbulent velocity of 0.2kms _1 were assumed. The data were 
convolved with a 6.7" beam, so that the disk fills the beam pre- 
cisely. For molecules containing 13 C and 18 the abundances of 
the main isotopic variant were used, scaled with factors of 1/45 
and 1/500, respectively. This assumption does not hold every- 
where in the disk due to isotope selective photodissociation, but 
the effects on integrated intensities over the entire disk are likely 
small. The resulting emission line profiles for O, CO, HCN and 
CN for models B1/BL1, B4 and BL4 are shown in Figure[9] 

It can be seen that the [O i] line does not change much when 
the mass of the disk decreases; the difference in intensity be- 
tween models Bl and B4 is only a factor of a few, while the mass 
changes with a factor of 1000. Model BL4 has a slightly higher 
[O i] intensity, since at the temperatures where O can be excited, 
less oxygen is in the form of CO and O2. Also, the [0 1] line has 
broader wings since it is excited mostly in the inner regions of 
the disk compared to model B 1 . 

The CO 7 = 6-5 line follows similar trends: it is opti- 
cally thick in all models, making it relatively insensitive to total 
disk mass. The line is significantly broader in model B4 than in 
models B 1 and BL4 because the CO column densities are more 
centrally concentrated in that model. 

The HCN 7 = 4-3 line shows a strong dependence on the 
disk mass. Model B4 has very low HCN abundances, giving this 
line an intensity in the mK range. The differences between mod- 
els B 1 and BL4 are less pronounced and reflect the slightly lower 



HCN abundances in BL4 with respect to Bl. The CN 7 = 3 - 2 
line, however, shows a strong dependence of the line shape on 
the dust/gas ratio. This line has a relatively low excitation tem- 
perature (~ 30 K), and the gas temperature in the region where it 
can be excited (i.e. where the density is equal to the critical den- 
sity of 10 7 cm 4 ) is too high in model BL4 to get an appreciable 
population in the upper level. In model B 1 the temperatures are 
lower due to the decreased C-ionization heating rate, and the 
7 = 3 level has a higher population. Model B4 only has a weak 
line here due to the low gas content of the disk. 

The integrated intensities are presented in Figure [10] It can 
be seen that model B1/BL1 generally produces strong emission 
lines, especially for the fine-structure transitions of O and the 
rotational lines of CO. These lines are the main coolants of the 
PDR surface of the disk; therefore they radiate all the energy that 
is delivered to the gas via the various heating processes. 

The [C 1] intensities show an initial increase in the B series 
even though the disk mass decreases by a factor 10, this is due to 
the high C abundance in model B2 compared to B 1 . In models 
B3 and B4 the intensity drops again, albeit slowly due to the 
optical thickness of the lines. In the BL series the [C 1] intensities 
only increase due to the ever increasing C abundances in this 
series. Similar effects also determine the intensities of the [Oi] 
line. 

The [Cn] fine-structure line is significantly weaker than 
those of [Oi] and [Ci], even though it is an important cooling 
line in PDRs. This line has a low critical density (10 3 cm" 3 ), 
which means is can only be an effective coolant at the immedi- 
ate surface of the disk. Deeper into the disk this line is quenched 
by collisions. Its intensity increases in the BL series because the 
diffuse upper layers are thicker there. 

The intensities of the CO rotational lines show a weak de- 
pendence on disk mass. Initially these lines are optically thick, 
but become optically thin going from model B3 to B4, causing 
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Fig. 10. Integrated intensities in a 6.7" beam for the B series (left column) and the BL series (right column). The top row gives the 
[C i] fine structure lines (7=1-0 solid, J = 2-1 dotted). The second row gives the fine structure lines of [C n] fine structure line 
(solid) and [O i] (7 = 1 - 2 dotted, 7 = 0-1 dashed). The third, fourth and fifth rows give the rotational lines of CO, 13 CO and C ls O 
(7 = 2-1 solid, 7 = 3-2 dotted, 7 = 6-5 dashed). The sixth row gives the rotational lines of HCO + (7=1-0 solid, 7 = 3-2 
dotted, 7 = 4-3 dashed). The seventh row gives the rotational lines of HCN (7=1-0 solid, 7 = 3-2 dotted, 7 = 4-3 dashed). 
The bottom row gives the rotational lines of CN (7=1-0 solid, 7 = 2-1 dotted, 7 = 3-2 dashed). 
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Fig. 11. Integrated intensities of HCN in the B series when all 
HCN is in the gas phase (solid line) and when HCN is removed 
from the gas phase for 7d us t < 70 K (dotted line). The BL series 
follows similar trends. 



a sudden drop in intensity. The lines do not change significantly 
in the BL series. 

The pure rotational lines of HCO + , HCN and CN follow the 
same trends as CO: they change only little in the first steps of 
the B series, but become optically thin going from model B2 to 
B3 due to their lower abundance compared to CO. These lines 
change very little in the BL series, with the exception of HCO + , 
which has a steadily decreasing column density in this series. 

3.5. Comparison with observations 

The rotational lines of CO, HCO+, CN and HCN of two Herbig 
Ae stars HP 1 63296 and MWC 480, have been observed by 
iThi et all ( 120041) with IRAM and JCMT The SEDs of these ob- 
jects indicate that they are group II sources, meaning that of the 
models presented here only steps 3 or 4 of each series need to 
be considered. Correcting for beam dilution, their intensities of 
CO, 13 CO and CN, as well as the CO/ 13 CO line ratio, fit well 
to the values found for model B2 (i.e. a 0.01 M Q disk with a 
dust/gas ratio of 0.01) with the HD 163296 data indicating a 
slightly higher gas mass than MWC 480. A good fit is also given 
by models BL3 and BL4 (i.e. a 0.1 M Q disk with a dust/gas ratio 
of 10~ 4 - 10~ 5 ), with the HD 163296 data indicating a slightly 
lower dust/gas ratio than MWC 480. The masses agree within 
a factor of a few with those derived from millimeter continuum 
observations assuming a dust/gas ratio of 0.01 . None o f the mod- 
els pre sented here can reproduce the HCO + lines by IThi et al.l 
(120041) . which have intensities of > 1 K km s 1 . This may be due 
to the low HCO + abundances found here compared to other disk 
chemistry models. Since these objects are group II sources, it is 
unlikely that model B2 describes their structure. Instead, models 
BL3 and BL4 provide a good fit to the observational data. 

The intensities predicted here for the HCN 4-3 line of mod- 
els B 2, BL3 and BL4 (~ 5Kkms _I ) exceed the upper lim- 
its bv IThi etal.1 d2004h of 0.8 and 0.3 Kkms 1 for HD 163296 
and MWC 480, respectively (corrected for beam dilution). This 
discrepancy may be explained by the fact that our models take 
freeze-out of molecules into account only by adjusting the global 
gas-phase elemental abundances, meaning that there is relatively 
little oxygen (since most of the dust is colder than the water sub- 



limation temperature of ~ 100 K) compared to carbon (since all 
dust is warmer than the CO sublimation temperature of ~ 20 
K). For nitrogen-bearing species this method is less accurate, 
however, since molecules that occur mostly in the disk's surface 
layers (like CN) will not freeze out due to the high dust tem- 
peratures there, while molecules that occur mostly in the deeper 
regions near the midplane (like HCN) are more likely to freeze 
out on the cold dust grains. To estimate the effect of freeze-out, 
the intensities were re-calculated by removing the HCN from re- 
gions where TVst < 70 K (consistent w ith the TPD experiments 
on H2S and C2H2 by Collings et al. 2004) . The results are given 
in FigureQj] It can be seen that freeze-out may account for more 
than an order of magnitude difference in intensity, and that the 
revised B2 model intensity is just consistent with the observed 
upper limits. The CN lines are less affected by freeze-out, since 
the peak in the vertical CN distribution occurs at larger heights 
than the HCN distribution, where the higher dust temperatures 
and lower densities will prevent freeze-out. 

The intensities presented here indicate that future instru- 
ments like ALMA will be able to detect and image the rota- 
tional lines of CO, 13 CO and CN in disks with masses as low as 
1O~ 4 M , with a typical sensitivity of 0.3Kkms _1 in 1 hr. The 
lines of C 18 0, HCO + and HCN can also be imaged in these low- 
mass disks, but will require long integration times. Massive disks 
that undergo dust growth and settling will be readily observable, 
since the molecular column densities remain high even when the 
disk becomes more and more optically thin in the continuum. 

There are several line ratios that can be used to ascertain 
whether a group II source (as determined by its SED) is the prod- 
uct of mass loss (putting it in the B series) or of dust growth 
and settling (putting it in the BL series). Some relevant ratios 
are plotted in Figure[[2] a low ratio (< 10) of CO/ 13 CO (< 10) 
or CO/C ls O would indicate an object with a low dust/gas ra- 
tio instead of a low overall mass; 12 C ie O is optically thick (thus 
making its lines relatively insensitive to the disk mass), while 
the isotopic variants of CO are optically thin, so they trace the 
CO abundance directly. 13 CO 3-2/HCO + 4-3 is another suitable 
tracer, since the HCO + intensities behave similarly in the B and 
BL series; on the other hand, the intensities of 13 CO drop off if 
the disk mass is decreased, while they remain high in the case 
of dust growth/settling. A high ratio (> 100) would indicate dust 
growth and settling, while a low ratio would indicate mass loss. 
The ratio of [0 1] 63 //m/146 fim directly probes the temperature 
structure of the upper layes, and is therefore suitable to distin- 
guish between the B and BL series. A high ratio (> 2) would 
indicate that a disk with a dust/gas ratio close to 0.01 ; a lower ra- 
tio would indicate dust growth or settling. The ratio of [Ci] 610 
jt/m/ 13 CO 3-2 can be used to distinguish group I sources (B1-B2 
and BL1-BL2) from group II sources (B3-B4 and BL3-BL4), 
where a high ratio (> 5) would indicate a group II source. This 
ratio is not very sensitive to the differences between the B and 
BL series, however. The CO 6-5/2-1 ratio was found to be insen- 
sitive to the differences between the B and BL series, since the 
gas temperatures stay high enough to excite the J = 6 even in 
models with a very low dust/gas ratio (e.g. model BL4). 



4. Conclusions 

The thermal balance, chemistry and emission line profiles were 
calculated for a series of disk structure models representing dif- 
ferent classes of protoplanetary disks around Herbig Ae stars 
with systematically decreasing disk masses or dust/gas ratios. 
The following can be concluded: 
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Fig. 12. Line ratios for the B series (solid lines) and the BL series (dotted lines). 



- For an interstellar dust/gas ratio (0.01), the gas temperatures 
at the flaring edge of the disk are very high, reaching tem- 
peratures of several thousand K. Directly underneath this hot 
surface layer is a warm (several hundreds of K) optically thin 
layer where the gas is heated directly by the stellar radiation. 
This warm layer drives an endothermic chemistry resulting 
in a relatively high CO abundance. In the midplane and in 
the outer disk the self-shadowing of the disk results in lower 
gas temperatures, usually close to the local dust temperature. 
This relatively cool self-shadowed part decreases in size with 
increasing disk mass. 

- The disk chemistry shows a strong dependence on disk mass. 
The low-mass disks have large optically thin atmospheres 
where fragile species like C2H and HCN are easily destroyed 
by UV radiation. More robust species like CO and CN still 
have relatively high abundances here. 

- Although the self-shielding of CO is stronger in the BL se- 
ries compared to the B series, the abundance of this molecule 
is lower. This is due to the reduced efficiency in CO forma- 
tion in the absence of PAHsand classical grains, and the high 
photodissociation rates of CO precursors like CH. 

- Dust settling has little effect on molecular abundances if the 
total gas mass is kept constant: there is only a small shift 
of the maximum abundance towards the midplane. The pri- 
mary effect of dust settling on the disk is to lower the gas 
temperature by decreasing the importance of the photoelec- 
tric heating process. 

- The CO rotational lines can be used to trace the total disk 
mass only for low-mass disks. In high-mass disks these lines 



are optically thick and not sensitive to the disk mass. The iso- 
topic variants of CO ( 13 CO and C ls O) are much more sensi- 
tive even at high disk masses; their intensities are also mildly 
sensitive to the dust/gas ratio. 

- Good tracers to distinguish between the B series and the BL 
series (i.e. whether a class II source is undergoing mass loss 
or dust growth and settling) are CO/ 13 CO and CO/C ls O line 
ratios, since the isotopic variants of CO are more sensitive 
to the total amount of CO in the disk than the main form. 
Also, the 13 CO 3-2/HCO+ 4-3 and the [Oi] 63 yum/146 yum 
ratios can provide information on whether a source belongs 
in the B or the BL series. [Ci]/ 13 CO can be used to distin- 
guish group I from group II sources, but is not sensitive to 
the differences between a disk undergoing mass loss or dust 
growth and settling. 

- In the disk models presented here, the dust temperatures are 
high enough to prevent the freeze-out of CO onto dust grains. 
For other molecules, like HCN, freeze-out remains possible, 
but has not been included explicitly in the current code. 
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